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Abstract 



In this paper we consider a class of the 2D integrable models. These models are higher spin 
XXZ chains with an extra condition of the commensurability between spin and anisotropy. 
The mathematics underlying this commensurability is provided by the quantum groups with 
deformation parameter being an A^th root of unity. Our discussion covers a range of topics 
including new integrable deformations, thermodynamics, conformal behaviour, S'-matrices and 
magnetization. The emerging picture strongly depends on the A^-parity. For the N even case at 
the commensurable point, ^'-matrices factorize into N = 2 supersymmetric Sine-Gordon matrix 
and an RSOS piece. The physics of an N odd case is rather different. Here, the supersymmetry 
does not manifest itself and the bootstrap hypothesis fails. Away from the commensurable 
point, we find an unusual magnetic behaviour. The magnetization of our chains depends on 
the sign of the external magnetic field. 



Introduction and Discussion 



Symmetry is the driving concept in particle physics. In Quantum Field Theory the particles 
are defined as finite dimensional irreps of the space-time and internal symmetry groups. In 
Statistical Mechanics the notion of symmetry has also played a very important role in the past, 
as a way of characterizing degrees of freedom and types of interaction. The recently introduced 
quantum groups are another stride in this direction, which deepens our understanding of sym- 
metry in systems with an infinite number of degrees of freedom. Quantum groups shed new 
light on the number of difficult problems which are, essentially, non perturbative in nature; 
determination of the particle spectrum, scattering matrices, correlation functions, to name just 
a few. 

An important family of integrable models is based on the quantum affine algebras Uq{G) 
[1]. These models are characterized by the quantum deformation parameter q and a finite 
dimensional irrep of G. The famous six- vertex model corresponds to G = Sl(2) and the 
fundamental spin 1/2 irrep. 0. The higher spin versions of Uq{Sl{2)) 0, can be constructed 
by fusion procedure 0. Generally, the spectrum of these models is expected to satisfy the 
bootstrap axioms |0: there is a fundamental particle such that all others can be interpreted 
as bound states. For the isotropic models of spin j (g = 1) the particle spectrum consists 
of a fundamental spin 1/2 particle which has (for j > 1/2) extra hidden RSOS spin [^], 0]. 
The corresponding physical S'-matrix factorizes into the product of the XXX S'-matrix and an 
RSOS piece. This picture remains essentially unchanged in the semiclassical weak anisotropic 
regime 0,0. The central extension for these nearly isotropic models is c = i.e. that of a 
SU{2)k WZW model with level k = 2j. Moreover, the excitations above the ground state can 



be understood in terms of the free bosonic and Z2j-parafermionic degrees of freedom |TT 



There is a connection between N = 2 integrable field theories in two dimensions and quan- 
tum affine algebras at roots of unity. This became clear in the study of the Sine-Gordon model 
in ref. |]12[ and more recently in |jl3|]. The N = 2 structure of the Sine-Gordon models appears 
at the special value = 87r3/2 which yields = 1 ( g = e*'^/^ , p = /3^/(87r — with the 
SUSY realized in a non local way W^, namely, the one defined by the non local action of the 
quantum group generators [|15|. The solitonic S'-matrix of the N = 2 theories factorizes into an 
= RSOS S'-matrix and the N = 2 S- matrix of the Sine-Gordon model at = 87r3/2 fTT 
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Observe that this factorization of the S'-matrices is similar to the one obtained for the physical 
S-matrices of higher spin XXX models using Bethe ansantz techniques. The only difference 
is that the XXX piece of the S-matrix is replaced in the case of the N = 2 models by the 
Sine-Gordon XX S-matrix (i.e. the one which corresponds to g"^ = 1). One could wonder at 
this point whether there exist spin chains whose S'-matrices factorize as those of the N = 2 
models. We shall give an example of these later. 

Commensurability for Uq{Sl{2))-models appears whenever the deformation parameter q is 
an A^th-root (g^ = 1), and the dimension of the irrep used in constructing the model is N' 
{N' = N {N/2) if odd (even)). For even the commensurable point determines the frontier 
of the weak anisotropic regime. The first result that we obtain comes from the explicit solution 
of the Bethe equations for the even case. From a technical point of view the novelty of the 
Bethe ansatz analysis is related to the fact that 2j -|- 1 is not a so called Takahashi number, an 



assumption that is done in previous studies of higher spin models 0, fl^, p6[. The elementary 
excitations S'-matrix has the same structure as that of A^ = 2 SUSY models, namely, a Sine- 
Gordon XX S'-matrix multiplied by an RSOS matrix associated with graph A2j+2- This fact 

ITT 

suggests that the commensurable spin j-chain with q = e^j+i provides a lattice model of A^ = 2 
integrable theories. 

There is another important reason to study the spin-anisotropy commensurable hamiltoni- 
ans, namely, the peculiarities of representation theory of quantum groups at root of unity |[T9| . 
At these values of q there appear new central elements, that is, new casimirs which form a 
closed Hopf subalgebra and whose continuous spectrum labels the irreps. Thus, one can change 
in a continuous way the irrep without changing its dimension. The models so obtained are the 
non-fusion descendants of the six- vertex model without classical q = 1 analog. The Chiral Potts 
model constructed from a class of irreps called cyclic, is one example of the latter [^. In this 
paper we shall limit our attention to a very special class of the non-classical irreps which we call 
nilpotent |2l|], |2^, |2^. These irreps have the highest and lowest weights. For the nilpotent 
deformations of the A^ even commensurable model we find a similar spectrum and the same 
central extension as for the commensurable point. The main difference is a replacement of the 
XX part of the physical S'-matrices by the nilpotent i?-matrix with A^ = 4. This nilpotent 
i?-matrix is related to an XX chain placed in the external magnetic field. However, the RSOS 
piece of the whole S'-matrix remains unchanged. The A^ = 4 nilpotent i?-matrix discussed 
above is, naturally, associated with the quantum algebra Uq{Gl{l, 1)), where q is fixed by the 
casimir which describes the nilpotent irreps. Since the GI{1, 1) algebra underlies N = 2 SUSY, 
we suggest that our nilpotent case is related to the g-deformed N = 2 SUSY. 

Another issue of interest appears while studying the behaviour of the nilpotent chains under 
the action of external magnetic field. The external magnetic field acts on all degrees of freedom, 
while the "nilpotent" magnetic field acts effectively only on the XX part. The net effect is 
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that the magnetization depends on the sign of the external magnetic field. The magnitude of 
this effect is proportional to the ratio of the XX degrees of freedom to the total number of 
degrees of freedom (i.e. it decreases for higher A^). 

We extend our analysis to commensurable spin chains with g-'^ = 1 and odd. The picture 
we obtain is drastically different from that of the N even case. The particle spectrum consists 
of the three types of particles whose energies do not satisfy bootstrap relations. The scattering 
shifts among these particles are the product of an RSOS piece associated with graph A2j^^.^.+2 
with jeff = j/2 (notice the renormalization of the spin) and another RSOS part associated 
with the "spin 1" model p = 2,r = 6 P3. The (p=2,r=6) RSOS model captures most of the 
aspects of our model: counting of states, central extensions and scattering shifts. However, 
there is a mismatch in energies, which reveals some failure of the bootstrap. For the odd 
case we find no trace of the N = 2 SUSY structure, but this is not surprising, since SUSY 
always requires an even number of degrees of freedom to be present. An attractive possibility 
will be to interpret the A^ odd case in terms of a Zs-parafermionic extension of A^ = 2 SUSY 

m. 



Let us summarize the main results we obtained: i) Integrable deformations of the spin chains 
away from the commensurable point. 

ii) Extension of the Bethe ansantz analysis to the cases where 2s + 1 is not a Takahashi number. 

iii) Computation of the S-matrices for the A^ even and odd cases, and connection with N = 2 
SUSY in the A^ even case. 

iv) New method of counting degrees of freedom in the vicinity of the ground state. 

v) Unusual Magnetic properties of the nilpotent models due to a non trivial interplay between 
an "internal" magnetic field associated with the nilpotent casimirs and the external magnetic 
field. 

The organization of the paper is as follows. In section 1 we present the basic tools for solving 
the descendent models including a brief description of the algebraic Bethe ansantz. Section 2 
contains some basics on the representation theory of quantum groups at roots of unity and a 
detailed description of the symmetries, i?-matrices and hamiltonians based on the nilpotent 
irreps. In section 3 we discuss the string hypothesis and solve the Bethe equations for the 
values of the nilpotent parameter, which yield the hermitian hamiltonians. In section 4 we 
study the A^ even case, computing the central extensions and the S'-matrices with N = 2 SUSY 
structure. In section 5 we extend this analysis to the N odd case. Finally, in section 6 we 
study the magnetic properties which depend on the direction of the external magnetic field. 
The paper ends with concluding remarks and two appendices where technical details are given. 
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1. Six- vertex Model Descendants 



1.1 Descendent Models 

Our starting point will be the quantum affine algebra Uq{Sl{2)) = Uq{A'i'^) [||]. The finite 
dimensional irreps of this algebra at level zero can be easily obtained from the finite dimensional 
irreps of the non affine quantum algebra Uq{Sl(2)) whose generators E, F and K satisfy the 
well known relation: 

EK = q^KE- FK = q-^KF; [E, F] = ^ ~ (1.1) 

q-q 

Denoting by irf an irrep of Uq{Sl{2)) of spin j of positive or negative spin parity we have: 

7lf{E)er = [r] Cr-l 

nf{F)er = ±[2j-r]e,+i (1.2) 
7rf{K)er = ±q^'-^'-er 

where r = 0, 1, . . . , 2j. This representation can be promoted to a finite dimensional irrep tt^^^ 
of Uq{Sl{2)) as follows: 

7rJ,(Eo) = e^vrf (F) vrJjEi) = e^vrf (E) 
7r^^„(Fo) = e-"7rf (E) vrj^jFi) = e~-Tif{F) (1.3) 
TrJjiTo) = vrf (ir-i) t^UK,) = 7rf{K) 

where u is called the spectral parameter. Given two irreps vr^j and 'n'j2,u2 the intertwiner 
matrix R^^'^^{ui — U2) is defined by: 

R'^'^^{ui - U2) TTj^^ui ® 7rj2,«2^(«) = ^ji,ui ® 7rj,„u2^'{a')R'^'^^{ui - U2) (1.4) 
The intertwiners (|1.4| ) satisfy the Yang-Baxter equation with spectral parameter: 

i?f/^(ui-W2) i?f3'^'^(«l-n3) I^^/\U2-U,) = RinU2-U,) i^fg'^'^ («i - Ug) i?^23''^ (^i -^2) (1.5) 
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In fact, the intertwiner i?-matrices, obtained through solving ( |1.4|) , are the trigonometric 
solutions to the Yang-Baxter equation ( |1.5|) and give rise to the integrable vertex models if one 
identifies the Boltzmann weights with the intertwiners R^^'^^ The simplest vertex model where 
the previous construction applies is nothing but the six-vertex model which corresponds to 
the choice ji = j2 = 1/2. Taking into account the fact that the spin 1/2 is a fundamental 
irrep of Uq{Sl{2)), one calls the {ji, j2)-w.odels descendants of the six-vertex model. Given 



the intertwiner W^'^^{u) one proceeds, in the spirit of ||2^, to construct the monodromy and 
transfer matrices: 

i) Monodromy matrix: 

Tn,n^u) = R^''^%u)a,LR'''''{u)a,L-l " ■ ■ R^''^' {u)a,l (1-6) 

where L is the number of sites of the chain and T^^'^^ is a (2ji + 1) x (2ji + 1) matrix which 
belongs to End{iS)^V^^) . 

ii) Transfer matrix: 

= Tr{T^''^^{u)) (1.7) 



From the Yang-Baxter equation ( |1.5| ), which is a consequence of the quasitriangularity of 



Uq{Sl{2)), one gets the commutativity of the transfer matrices: 

t^''^{ui),t^''^{u2)] =0 (1.8) 

which implies, in particular, that the transfer matrices of the (j, j)-model can be simultaneously 
diagonalized together with the one of the (1/2, j)-model. The main advantage of this result is 
that the algebraic Bethe ansantz can be easily extended to the (j, j)-models. 



1.2 Algebraic Bethe Ansantz for Higher Spin 

Let us, first, consider the case of a quantum deformation parameter q not equal to the root 
of unity. We shall briefly summarize the algebraic Bethe ansantz for a vertex model of spin j. 
The basic tools that one needs are : 
i) Reference state: 



/ 1 \ 

^0 



;i.9) 



^ Writing u as hu and q as e'' and taking the limit h ^ one gets the rational solutions depending on u. 
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ii) Spin-wave creation operator: 



Bj{u) 



u 



1,0 



;i.io) 



Recall that T^^'^'^{u) is a 2x2 matrix. The action of the Bj{u) operator on the reference state 
(|1.9|) reduces the total spin of this state by one and creates a spin wave whose quasimomenta 
depends on the spectral parameter u. 

iii) Spin-wave annihilation operator: 



0,1 



;i.ii) 



Since the spin j irrep is a highest weight vector irrep it follows: 



Cj{u)\Qj>=0 (1.12) 

Thanks to (|1.8|) the candidates for the eigenvectors of the transfer matrix t^/^'-^ are of the 
form UT=iBj{ui): 

m m 

iMu) + D,iu)) l[BM)\^^ >=A,iu,u.)l[B,{umj > (1-13) 

i=l i=l 

where 

Mu) = {T'/'Hu))o,o , D,{u) = {T'/''^{u)),,, (1.14) 

The conditions on the spin-wave rapidities Ui imposed by ( |1.13D can be easily obtained using 
the commutation relations between the creation operators Bj{u) and Aj{u) and Dj{u). These 
commutation relations are encoded in the famous RTT=TTR equation which follows from the 
definition of the monodromy operators and the Yang-Baxter relation ( [1.5| ): 

^iij rj.ni j^hj ^ j^hj j^jij ^in ^^_^5) 

where Ti = T ® 1, Ts = 1 ® Tg. 

Taking ji = js = 1/2, j2 = j, we get the desired commutation relations. Indeed, equation 
(|1.15|) defines a so called Yang-Baxter algebra whose "structure constants" are determined by 
the six- vertex i?-matrix i?^/^'^/^. Since the irrep vrj is a h.w.v. representation, one deduces that 
Aj{u) and Dj{u) are (2j+l) x (2j+l) diagonal matrices. Defining aj{u) and bj{u) as: 



Mu)\nj >= iajiu)f\nj > , B,{u)\n, >= {b,{u)f\n, > (i.ie) 

the Bethe equations for the rapidities Ui read: 
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/ aj{ui) \ ^ H a{uk -ui)h{ui - Uk) 
\bj{ui) J f^Jll^^ a{ui - Uk)b{uk - ui) 

where a{u) = ai/2{u) and b{u) = 61/2 (m). Notice that only the Lh.s. of equation (|1.17| ) depends 



on j, the r.h.s., corresponding to the six- vertex model, is the same for the whole family of 
descendants. It will be convenient for later use to introduce the following change of variables: 
u ^ ^{u — i). In terms of these new variables, the Bethe equations become: 

/ sinh^jui + 2j _ |1 sinh^{ui - Uk + 2i) 
\sinh^{ui - 2j i) J f^J^^^,^^ sinh^{ui - Uk - 2i) 

The hamiltonian of the spin chain is obtained through the logarithm derivative of the transfer 
matrix: 

H,{-f)=tI^logt"{u)\u=o (1.19) 
where / is an overall coupling constant. 

For j = 1/2 one gets the well known XXZ-chain with anisotropy A = ^±2 — = cos'y. In the 
limit, when the parameter 7 goes to zero, one recovers from ( |1.19| ) the higher spin isotropic 



models of the first two refs. in An example of the spin 1 isotropic hamiltonian is: 

L 

i?,=i(7 = 0) = ^ ■ - {S,S,^if (1.20) 

i=l 

Within the family of anisotropic higher spin integrable models we shall pay particular at- 
tention to the cases where the anisotropy 7 and the spin j are commensurable in the sense: 



= 1, 2j + l = N, N odd 

= 1, 2j + 1 = N/2, N even (1.21) 

where q = e^"' . At these points, the higher spin XX Z models admit continuous integrable 
deformations parametrized by the eigenvalues of the central elements of Uq{Sl{2)) at = 1. 
This will be the subject of the next section. 
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2. New integrable Deformations for 

N 1 
q^^ = 1 



2.1 Mathematical Preliminaries: Representation The- 
ory of (7,(51(2)) 

The representation theory of quantum groups at roots of unity has been worked out in full 
generality in refs. |TB[, |]TP|. Here, we will only need the results concerning Uq{Sl{2)), which 
are summarized below. 

For = 1 the Hopf algebra Uq{Sl{2)) contains a central Hopf subalgebra Zq generated by 
the elements E^\F^' and where N' = N {N/2) if is odd (even). The Schur lemma 
implies that the finite dimensional irreps tt are characterized by the eigenvalues of the central 
elements: 

7r(a) = Ua)l , a E Zq (2.1) 

The finite dimensional irreps can then be divided, according to the eigenvalues of the central 
elements, into the following four types: 



irrep tt 




UE^ 




dimension 


hwv / Iwv 


cyclic 




Vt^O 




A^' 


no/no 


semicyclic 







2 ^ ±1 


A^' 


yes/ no 


nilpotent 








generic 


A^' 


yes/yes 


classical 








±1 


< A^' 


yes/yes 



Table 1. 



The last two columns refiect the highest and the lowest weight properties of the irreps. The 
classical irreps are often called regular or restricted. We prefer to use the name " classical" , in 
a sense, that they originate from the continuous deformation of the finite dimensional irreps of 
the classical algebra. In what follows, we shall concentrate on the nilpotent representations. 
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The method of section 1 for constructing the integrable vertex models associated to irreps of 
quantum groups, can, automatically, be applied to the nilpotent ones. These representations, 
that we shall denote by ttx, are given in a basis {cr}^ ~^ by: 

nx{E)er = dr-i{X)er-i 

7ix{F)er = d,_i(A)e,+i (2.2) 
TX\{K)er = Xq~'^'^er 

where 



dl{X) = [r + l] \'-'-\''' (2.3) 



q-q 



The eigenvalue of the central element K^' is A^'. The representation ( |2.3| ) has been chosen 
because the raising operator is the transpose of the lowering one, i.e. (nx^E))^ = TTxiF). 
The complex parameter X E C labels the nilpotent irreps. It is convenient to introduce a 
"generalized" spin s by the relation 0: 

A = q^', 2s = 2s modN (2.4) 



For generic values of A, the representation (|2.3| ) is irreducible. There are, however, a finite 



collection of A's for which the representation is reducible and fully reducible into the direct sum 
of two classical irreps. This happens for the following values of A: 



AL^)=±g^ m = 0,l,---,iV'-2 (2.5) 



for which 



^aL±' = ^™/2 ^(^'-m-2)/2 (2-6) 

We will call the values (|2.5| ) orbifold points. Note that they are fixed points under the action 
of the coadjoint group defined in . 

Next we want to study certain automorphisms acting on the space of the nilpotent irreps. 
Let us introduce two discrete transformations called spin-parity (P) and charge conjugation or 
spin-flip (C) as follows: 

A^ = -A , X^ = q-^X-^ (2.7) 



Throughout this paper we use j and s to designate "regular" and nilpotent spins, respectively. 
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From (|2.3| ) one derives the relation between the corresponding irreps: 



TT^p (E) = iTTx (E) , Ti^c (E) = Cnx {F)C 

n^p{F) = tn^{F) , n^ciF) = Ctix{E)C (2.8) 
7T^p{K) = -7i^{K) , n^c{K) = Cn^{K-')C 

The charge conjugation matrix C is defined by: 

C (Cr) = CN'-r-l (2.9) 

The proof of (0|) is based on the relation dr{X'") = (iAr/_r-2(A). A representation tta is 



called hermitian if lEBI: 



nxiEY = vxnxiF) , 7ix{K)^ = TTxiR-') (2.10) 

where vx = ±1 is a spin-parity of the irrep. Equation ( |2.10|) implies that A is a phase (i.e. s is 
real) and that the spin parity depends on the spin s as: 

/ sin2js\ 

Vs = vx = sign — (2.11) 

\^ stn'f J 

In general, equation ( p.lO|) implies certain restrictions on the allowed values of A : 

vxdl{\)>0 for r = 0,1,...,N' -2 (2.12) 

or equivalently 

Vs sin-^k sin-^{2s + 1 - A;) > /or A; = 1, . . . , iV' - 1 (2.13) 

Hence, for a hermitian irrep the sign of d^{X) is independent of r. Our convention for taking 
the square root is: 

4(A) = w ''>;'+", r.!;',-Lxi/2 (2-14) 



A hermitian representation tta of Uq{Sl{2)), when lifted to a representation nx^u of Uq{Sl{2)), 
satisfies the following relations: 

[7rA,„(^o)]^ = ^^A 7rx,u{Ei) 

[t^xAFo)]^ = vrA,„(Fi) (2.15) 
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where we have assumed that u is reaL 

These relations are preserved by the comultiphcation. The behaviour of the spin-parity under 
the discrete symmetries ( |2.7| ) is: 



vxp = -vx , vxc = vx (2.16) 



From ( |2.13| ) we deduce that the allowed values of the spin s belong to open intervals whose 



boundaries are the orbifold points ( |2.6| ). From now on, for the sake of simplicity, we deal only 
with cases where q = e^'^''-/^ , The hermiticity intervals for the cases of even and odd are 
given in figure 1. We observe that the P-operation interchanges the two spin-parity intervals 
(i.e. 2s ^ 2s + N/2), while the C-operation maps each spin-parity interval onto itself leaving 
invariant the middle point (i.e. 2s — — 2s — 2). In fact, the middle points of these intervals 
correspond to classical higher spin irreps whose spin j is commensurable with the denominator 



of the anisotropy, as in equations ( |1.21| ). As can be seen from table 1, these irreps are both 
classical and nilpotent. The boundaries of the hermiticity regions for q = e^'^*/^ are the following 
orbifold points: 



N even : X= ±l,±g2 

N — 3 N — 1 

Nodd: X= ±q—,±q— (2.17) 

which together with ( |2.6| ) impliy that at the orbifold points the hermitian nilpotent irreps break 
according to the pattern 1 © ^ for A^ even and © for A^ odd. In appendix A we study, 
in full generality, the hermiticity regions associated with generic g-primitive root of unity. 



2.2 q-Chains Based on Nilpotent Irreps 



In this section, we want to find the intertwiner i?-matrices associated to the tensor product 
'^\i,ui ® '^X2,U2- These _R- matrices can be found, using a recursive method for solving the 



intertwiner equations ( |1.4| ) in the case of nilpotent irreps. Taking into account the normalization 
condition 



pOO 
-^00 



(2.18) 



we have 



U-lV''~\e-X,X2q' 



■X 
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ri r-2 
Zi=0 «2=0 



ri 




^2 


. ^1 . 




. ^2 . 



■{q-q- 



^■1+^2 — 1 ri+h — l r2 — l ri+r2— Z— 1 

X n ^.(^i) n ^.(^i) n ^.(^2) n ^.(^2) 

j=ri j=li+l2 j=r2-l2 j=r2-l2 

r2-i2-l h-l 

X Al^A^i-'i n (e"A2g-^' - e-"Aig^) [] (e"Aig-^+''^-'2 - e-"A2g^+'^-''^) (2.19) 
i=o j=o 

This is another trigonometric solution to the Yang-Baxter equation. Strictly speaking, so- 
lution (|2.19| ) comes with a factor g'(''i+''2-0-»'i»'2 which can be eliminated under an appropriated 
rescaling of the basis. Next, we list some useful properties of the i?-matrix (|2.19|) : 

i) Normalization 



RxAu = 0) = P 

where P is the permutation operator P{v i^w) = w0v,v,wE C^' . 
ii) Unitarity 



(2.20) 



iii) Parity 



iv) Spin-Parity 



R\^,X2{u)Rx-,^X2{-u) = 1 



P Rx^,X2{u)P = i?A2,Ai(M) 



(2.21) 



(2.22) 



Rxp,x-iu) = (1 ® e^'^^)i?,„,,(«)(e*-^ © 1) 
where A/" is the number operator: 



(2.23) 



AT. 



v) Charge- Conjugation 



1 



Rxc,xc{u) = -^i^v'-l., . 

-Kjv'_i,Ar'_il^l, ^2, U) 



vi) Hermiticity 



(2.24) 



iC0C)Rx„x2iu)iC®C) (2.25) 



^La2(") = ^Ai,A2(-m) 



(2.26) 
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This last equation holds true when both irreps Ai and A2 are hermitian and have the same 
spin parity. It reflects the fact that the dagger operation is compatible with the comultiplication 
of Uq{Sl{2)), namely: 



®vrA2,„2(A(Fo))]^ = Vs, 7rA„„, ®7rA2,„2(A(Fi)) (2.27) 



Using the i?-matrix ( p.l9|) we can deflne an integrable vertex model whose integrability 
condition Hi. 81) reads: 



[t'^'\m),t'^'\u2)] = (2.28) 
In this case, the Bethe equations are a generalization of ( |1.18| ) to continuous values of the 



spm: 

/ sinh^juj + 2si) \^ _ ^ sinh^{uj -Uk + 2i) 
\sinh^{uj - 2si) j ^=1,1;^^ ■5«'^/i^(uj - - 2z) 

The spin-chain hamiltonian ( |1.19| ) will, of course, depend on the spin s: 

Hx = Hsii) = ^I-^log t^'\n)|„=o (2.30) 

where / is an overall coupling constant. 

Before giving the general form of these hamiltonians, it is worth to exhibit their properties 
which follow from those of the i?-matrix: 

i) Hermiticity: If vr^ is hermitian then 

Hi = Hx (2.31) 

ii) Parity 

PlHxPl = Hx (2.32) 

where Pl is the operator that reverses the order of the sites along the chain. 

iii) Spin-Parity 

Hxp = nHx^^ (2.33) 

where n = nj=i e^^^'^-'^ . 

iv) Charge-Conjugation 

Hxc = VHxV -zL^R^:zl^'rl{X,X;u)\u=o (2.34) 
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where V = Il}=iCj. 

Since Q and V are unitary operators, we conclude that the hamiltonians associated to X,Xp 
and Ac describe the same physics. Hence, we shall restrict ourselves to a fundamental region 
of the P and C operations in the A space. For any value of q we can choose the half intervals 
going from the classical point up to the orbifold point within a hermiticity region. In the cases 
where q = e^'^*/^, these intervals are given by: 



N even : j < s < j + 1/2 
N odd: j <s<j + 1/4 (2.35) 

Thus, we have a one parameter family of Hamiltonians Hs{'~f = 2ti/N) on the interval ( p.35|) , 
which we can compare with the higher spin XX Z hamiltonians Hj{'y) of section 1. The results 
of this comparison are graphically presented in figure 2. Along the line I, one continuously 
modifies the anisotropy 7 keeping the spin j fixed, while along the line II, one varies the "spin" . 
These two lines meet at a single point A which corresponds exactly to the commensurable 



higher spin cases of ( 1.21 



A peculiar feature of the nilpotent models (i.e. line II) is the existence of another local 
conserved quantity obtained by taking the derivative of the transfer matrix with respect to Ai 



Qx = 2X^^hgt'^'\u = 0)\ 



Ai=A 



The counterparts of (|2.21| - p.26|) for Qx are: 
i) Hermiticity 



^2.36) 



ii) Parity 



iii) Spin-Parity 



iv) Charge-Conjugation 



PlQxPl = -Q> 



Qxp = -nQxn^ 



(2.37) 



(2.38) 



^2.39) 



Qxc = -VQxV (2.40) 

Next, we shall consider a few examples in order to get some insight into the properties of 
the models we intend to study. 
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Case N=4 



The point A of figure 2 has 7 = 7r/2 and j = 1/2 and corresponds simply to the XX modeh 



(2.41) 



i=l 



which is equivalent, after a Jordan- Wigner transform, to a free fermion model. Line I of figure 
2 corresponds to the XX Z models for different values of the anisotropy 7 . The hamiltonians 
on line II are: 



A - A 



-E( 



(A + A ^) z I \\ 



(2.42) 



i=l 



and can be recognized as those of the XX models with an external magnetic field given by 



A+A- 



cosTis. The hermiticity region of (|2.42| ) is 2s G (0, 2) whose boundaries correspond to 



the critical values of the magnetic field. The operator Q is given by: 



A- A 



■t=i 



X Y 



o, a. 



^y \ 



(2.43) 



This is a hopping hamiltonian which is parity odd in contrast with ( |2.42|) which is parity 
even. Combining these two commuting operators we can construct a hopping hamiltonian with 
a complex hopping parameter. The hamiltonian ( p. 42 ) generalizes the free fermion model in 



a quantum group symmetric way (the corresponding i?-matrix was first introduced |2^ from 
somewhat different point of view). In fact, the nilpotent i?-matrix coincides, in this case, with 



the i?-matrix of GZ(1, 1)^ (g = A) in the fundamental representation P9 |. 



Case N=3 

The hamiltonians on line I were constructed by Fateev and Zamolodchikov and are associated 
with the spin 1 XX Z models [Q: 

Hfz{q) = E "^fc ^k+l + ^k+l H 7) ^k^k+l 

k=l ^ 



2 2 

( cX qX j_ ny qy \2 Q l qZ qZ \2 



(2.44) 



+ (1 — q — q ^) {S^S^_^i + SX SX+i)S^ S^j^i+ 
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where 3^,3"^ and 3^ are the classical spin 1 matrices. At the isotropic point 7 = 0, ( 2.44|) 
becomes the hamiltonian (|1.20|) . 

The later model is gapless and corresponds to a WZW model 3U{2) with level k = 2j = 2 
which has central extension c = 3/2 0. The hamiltonians on line II can be written as: 



qY qY ^ 



" qZ qZ 



k=l 



i^k ^k+1 + ^k+lY + l^i^k ^k+l) 



■<Y cY 



+ ( - h dodij (3^ 3^^i + 3]^ 3]^^i)3^ 3i 



z 

k+r 



(2.45) 



\q 



2(g-g 



i^k ^k+1 + ^k+l)i^k + 'S'fc+l) 



2{q-q' 

which coincides with ( p.44| ) at point A , i.e. A = up to an overall factor. From this result 
we see that the hamiltonian Hpz{q) at = 1 possesses an additional symmetry given by the 
operator Qx with X = q^. 

In the = 4 case both hamiltonians Hx and Qx were seen to describe the hopping of a free 
fermion in the presence of a magnetic field. In the = 3 case one has three degrees of freedom 
per site. This fact can be accommodated in a model of two fermions A and B subject to the 
constraint of no double occupancy. Upon the identification cq ^ A,ei ^ 0, 62 — ^ B, where 
denotes an empty site, we deduce that the hamiltonian (|2.45|) at A = describes the following 
processes in the terminology of the ref. 



30 



1) Diffusion 



A + 0^0 + A, 5 + 0^0 + 5, A + B^B + A 
dS + A^A + Hs, dS + B^B + dS, B + A^A + B 



^Alh - - ^B(h - - -^AB - -J^BA - l2-4Dj 



where 7 = 2it/3. 
2) Annihilation 



stri'y 



A+B ^ 0+0 
B+A ^ 0+0 

hZ= H^b'a =-- ^ (2.47) 
stwy 
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3) Creation 



+ ^ A + B 
0+0 ^ B+A 
Hm = ^00^= -— ^ (2.48) 

The hamiltonian has in addition a chemical potential term /i 1 + 1 ® /U with : 

fj.A = fJ^B = TT- — At0 = (2.49) 

We obtained a "chemical interpretation" along the lines of ref. of the hamiltonian 

H^=,2 = Hpziq' = !)• 

Another interesting hamiltonian emerges when one reaches an orbifold point. As can be 
seen from ( |2.45| ), Hx has a simple pole at an orbifold point Xq. Then, it makes sense to define 



an orbifold hamiltonian H^^^ as a residue of H\ at this point, namely: 

N'-2 

= n dl{X)H, (2.50) 
Choosing Aq = 1 we get the following "orbifold" hamiltonian: 

HiZ? = E + rrrl^, + p,^, + p^^, - - a%,) (2.51) 

where T,p,a are different embeddings of Pauli matrices acting in C^: 

T- = E^\ p- = (2.52) 



and E^^ is the 3x3 matrix with 1 in the i — j position and zeros elsewhere. Eqn. ( p.51|) has a 



classical SU{2) x f/(l) invariance generated by J2j <^fi Yjj ■, Iii ref. this hamiltonian 



was given a chemical interpretation in terms of diffusion processes only. Another possible 
interpretation of ( |2.51| ) is a Hubbard model with infinite Coulomb repulsion or, equivalently, 
a t — J model with zero spin-spin coupling. It is known that these models have an infinitely 
degenerated ground state when L oo. This degeneracy is broken in the t — J model by the 
spin-spin coupling while, in our case, it is broken by the values of A different from 1. We will 
further discuss = 3 case in section 5. 
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Generic N Case 



Taking into account the complicated structure of the hamiltonian (p.45|) for the = 3 case, it 
would seem hopeless to find an explicit expression of the hamiltonian for the generic values of 
A^. Fortunately, this is not so and, in fact, we get the following expression: 



L N'~l 



1 



S+(A)ST^i(A)) + (ST(A)S+ (A)) 1 +;x,(A) +/i,+i(A) (2.53) 



where the matrices S=''(A) are given by: 

S+(A)e, = 
S-(A) 

and 



q' - q- 



1/2 



g'^-iA-1 - g-^'+iA 



Cr-l 



(S+(A) 



(2.54) 



/i(A) Co 

/i(A) Cr 



:^\-^qi + \q-i 

-1' / -; — — ■ ; ■ 

P,\-'q'-Xq-^ 



[2.55) 



In deriving (|2.53|) , we have assumed that A has a negative parity. 

The hamiltonian (|2.53|) exhibits strong formal similarities with the one of the Chiral Potts 



models 120 



L N-l 



j=l n=l 

where the couphng constants a„ and are given by: 



:2.56) 



e N 



oik 



sm 



' -irk \ ' 
, AT ) 



cosd) = k'cosc 



(2.57) 



and X and Z satisfy = = 1, XZ = e^'^'/^ZX. A possible basis for X and Z 



is: 



XCr 
ZCr 



Cr-i modN 



(2.58) 



To bring closer the relationship between ( |2.53|) and (|2.56|) we will consider the case X = q ^ 
and 7 = ir/N, namely: 
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1 




2n. 



(2.59) 

where = S^(A = q^^) are an A^-dimensional version of the Pauh matrices (indeed, 
= 0"^ for N = 2). The origin of the formal similarity between ( p.59| ) and the Chiral Potts 
hamiltonians is due to the chain of reductions cyclic — > semicyclic —>■ nilpotent irreps and the 
fact that the chiral Potts model is built up from the cyclic irreps. Indeed, if we choose = 
and 4> 7r/2 in ( p.56|) we get the same Z terms up to a divergent piece proportional to the 



number operator. Still, one has to explain why the cyclic operator X reduces to the nilpotent 
operator Although there are formal similarities between the Chiral Potts model and the 
one discussed here, we believe that these two models are not equivalent. The main reason is 
that Chiral Potts preserves a Z^ symmetry which is broken in our case. 

The discussion above suggests that the models we constructed belong to a more general 
class of integrable theories, whose hamiltonians have the generic structure: 

L N'-l 

HCN = -T.T. [tn {^t^J+lT + tn + (2.60) 

j=l n=l 

This kind of hamiltonians can be obtained through the extension of the Hopf algebra 
Uq{Sl{2)) by a central element z with non trivial comultiplications. In this way, the parity 
is broken and one gets chiral hamiltonians. 

To end this section, we want to analyze what happens when we reach an orbifold point of 
the hamiltonian (|2.53| ). If A approaches ±g™, we see that has the block diagonal form 



Hj'/j or HjYi, where I and II denotes the two sectors in which the nilpotent irrep breaks at 
the orbifold point. This means that the only processes taking place are diffusion ones between 
states which belong to different sectors I + II ^ II + I. 
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3. Thermodynamics 



3.1 Bethe Equations, The String Hypothesis 



In this section we will consider in some details the Bethe equations ( |2.29| ) with generic spin 
s in the hermiticity region ( p. 35 ). In the thermodynamic limit L oo, we will assume the 
string hypothesis |T^. This means that the solutions to ( p.29[ ) are organized into the strings of 
roots sharing a common real value. These strings are characterized by length n and parity f„ 
as follows: 

W , .-u , 1 _ o; , M . / = i,...,n (3.1) 



u 



u + i[n + 1 - 21 + —{1 - VsVn)] 
27 



where u represents common the real part and Vg is the spin parity introduced in (|2.11 ). Strings 
of roots can be intuitively interpreted as bound states of the spin waves. From the Bethe 
equations ( p.29| ) we derive the consistency conditions for the allowed strings ( |3.1| ): 



^ g/if(^"^ + 2zg) 
j=i sh^{uf'^ — 2is) 



> 1 ; k 



n — 1 



(3.2) 



It is not hard to prove that the consistency condition ( |3.2| ) for the Takahashi strings ( |3.1| ) is 
equivalent, for asymptotically large rapidities, to the hermiticity condition ( 2.13 ) for irreps of 
dimension n and spin parity f„. In fact, for m 00 we get from 



sm(27s) .... , , . 

Vg Vn sm(7n) sm7(n — k) > ; 

sin 7 



k 



,n — 1 



(3.3) 



Using ( |2.11| ) we obtain the hermiticity condition ( p.l3| ) for irreps of dimension n. Notice that 
the condition (|3.2| ), for m 00, fixes the allowed strings in a form completely independent of 
the value of the generalized spin s. To feel s we need to consider ( |3.2| ) for generic values of the 
rapidities. We now show that if s is in the hermiticity region ( |2.35|) then the asymptotically 
allowed strings ( p.3|) will satisfy ( p.2|) for arbitrary values of rapidities. To this end we rewrite 
fO) as: 



^ z - cosjy + 2-f{s - j)) ^ ^ 
z - cos{y - 27(5 + j)) 



(3.4) 
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where z = ch'yu; y = 7(n + 1 + Po^~~f^) and 1 < A; < (n — 1). It is straightforward to prove 
the following identity: 

-pj- 2; - cos(?/ + 27(5 - j)) ^ jl 2vnVsSm'y{2s - I) smj{n - k + I) sm'y{k - I) 

/ii 2;-cos(y-27(s + j)) ~ ^ i^o i=o sin7(/ + 1) z - VnVsCos-f{n + 1 - 2k + 21 - 2s) 

(3.5) 

Taking into account the hermiticity condition ( |2.13| ) and the Takahashi condition ( |3.3| ) one 
immediately concludes that l.h.s. of eqn. ( p.4|) is, indeed, greater than one. 

Making use of "Takahashi zone" terminology (also see appendix A) we have, for allowed 
strings {nj,Vj), the following: 



even : 



— zone rij = j, f„^. =+1 l<J<z^— 1 

1 — zone riy = l,v^ = —1 j = i' 



(3.6) 



N odd: 



— zone rij = j, f j = +1 1 < J < — 1 

Tiy = l^Vu = -I i = V 

n^+i = z/ + 1, v^+i = +1 j = z/ + 1 

2 - zone ny+2 = v-, = +1 j = v + 2 



1 — zone 



where u = ^i^^Y^) for A^-even (odd). 

In the thermodynamic limit, equations (|2.29| ) become: 

a,{u) = i-iy^ip, + pf) + Y.Tjk * Pk (3.7) 

k 

where Pjip^j^'') is density of j-strings (j-holes) and (— l)*"^ = sign(aj). The Fourier images of 
the functions which appear in the equation above are given by: 



^jk = g{w;\nj -nk\;Vn^VnJ + g{w;\nj + nk\,Vn^VnJ + 

min{nj ,nf^)~l 

+2(1 - Si^rnin{n,,nk)) di^'i l^j - + 21] f„^.W„J (3. 

1=1 



ttj = g{w] 2s + 1 — rij + 21, VsV. 



1=0 



sh2pow{{^ + ^)) 

g(w,n;v) = ; {(x)) is Dedekind function 

shpQW 

The "bare" energy of strings of the type j is given by: 
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An 

E„^{u) = -—la^iu) (3.9) 
where I is an overall coupling constant which is assumed to be positive throughout this paper. 



Clearly, the equation (|3.7| ) makes sense only if the sign of aj{u) (or Enj{u)) does not change 
when u varies from zero to infinity. It is a bit surprising that no proof of this fact has been 
given in the literature. The following is the proof that, indeed, there is no sign crossing for 
aj{u) as long as rij is an allowed Takahashi string length and s belongs to the hermiticity region 
(|2.35|) . We represent energy En in the following form: 

2IVnVs sin(27g) sin(7n) 
sm7[z — VnVs cos7(2s + n — 1)\ 
-LllX!n 2t;„t;,sin7(2g-/)sin7(n-Z)sin7/ ^ 



m- 



iJl sin7(/ + l)[z — VnVs cos7(2s + n — 2/ — 1)] 



Once again, making use of formulas (|2.13| ) and (|3.3|) , one can infer that the expression inside of 



the figure brackets in the equation above is always positive and, therefore, the sign of En is: 

sign(En) = — sz5f?T,[f„fs sin(27s) sin(7n) sin7] = —(—1)''^ = —(—1)* (3-11) 

with i being a label of the "Takahashi zone" where rij string lives. 

Before we move on, let us make the following comment. Irregular spin s, in some sense, 
plays the role of the "internal" magnetic field, producing, for instance, the nonvanishing ground 



state magnetization [^. However, this interpretation cannot be taken literary because other 
attributes of magnetic field are absent. For example, none of the string energies change their 
signs as opposite to the case where system is exposed to the external magnetic field. 
Renormalized or "dressed" energy of the string is defined as: 

exp [/3e,{u)] = = r/, ; P = - (3.12) 

Pj{U) 1 

with T being temperature. 

Equations for the "dressed" energies can be obtained by minimizing the free energy F 



/ + 00 
du\aj{u)\ln{l + rij^) (3.13) 
3 



F ^ /•+00 



in the presence of the Bethe ansatz constraints ( p.5| ) as was first proposed in ref. |31|. The 
resulting Dyson like equation is: 

Inr^, = -4po//3a, + ^(-l)^^T,-, * ln(l + r^^^) (3.14) 

k 
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Before giving the explicit solution to ( 3.14 ) let us briefly comment on its physical meaning. 
The positive energy strings (and the holes in distribution of negative energy strings) can be 
interpreted as particle like excitations. The negative energy strings will determine the ground 
state defined as a state where all negative energy modes are filled. Finally, the zero energy 
strings reflect the symmetries of the model and will be used to define the internal quantum 
numbers of all elementary excitations. 

To solve equation ( p.l4| ), it is convenient first to invert the matrix T (see appendix B). Then, 
this equation can be converted to the following form: 



Case N even 
In rij 

In r]^-i+f 

Case N odd 
In rij 

In rj^ 

In r]^+i+f 



si * [ln(l + r7,+i)(l + r],^i) + 5,-,„2 ln(l + Vu^)] ; 1 < J < - 2 

Ipo sin TTx 



[-lY[s,*\nil + r],_2) + 



T ch^ + (-l)^cos7rx^ 



f = OA 



(3.15) 



Si * [ln(l + r]j+i){l + r]j^i)] + 5j>_i[si * ln(l + r]j) + 



+ {s2-si)*\n{l + r]^) 



4po/ 
T 



/(«)] 



1 < J < «^ - 1 



S2 * ln(l + r^^+i) (1 + vM) i^+Vu-i 



1)^^^ -I- ln(l I ry ) cos27ra: f = 1 

T c/iTTu + (— 1)'' sin 27rx ' ' 



where si(u) 



7 seek 



2 ' 



/ = 1, 2; ?7o = 0; X = s - po 



3-Vs . 
4 ' 



The equations above can be easily solved in two cases: u ^ co (or /? ^ 0) and -u ~ |T7 . 
These cases will be important later for the central charge calculations. Results are summarized 
in table 2 below: 



even 


odd 


i + ^i = (j + i)'; i<j<i^-2 


i+% = (j + i)'; 1<J<'^-1 
l + ^. = ^'+i = (l + ^^')' 

= !-(! + jy)-' 


= oo 


, , = sin" ijXaO'+l) i ^ ■ ^ i 



Table 2. 

In table 2 above, = ri{oo) and = 1]{0). 

It is relatively simple to convince oneself, using equations ( p.l5[ ), that all strings are non- 
negative except z/ — 1 string for even and i/ + 2 for N odd. Making use of this fact, it is now 
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easy to solve equations (|3.7|) and equation ( p. 15 ) in the limit of zero temperature to obtain the 
following results presented in tables 3 and 4. 





Ground state strings 
(p^'^) = 0) 


Positive energy strings 
(P = 0) 


Zero energy strings 

(P = P^') = 0) 


even 


z/- 1 


V 


the rest 


odd 


I/ + 2 




the rest 


Table 3. 



Negative Energy 



Positive Energy 

sW(2s+2-po^^^) 



even 



4poJ 



-pv-\ 



shiw 



4po/ 



PI 



(h) 



sh2w 



odd 



<;i/+2(w) 

4po/ 



~ _ ry Chfshw{2s-po':^) 

-pv+2 - ^ ^f;^ 



Apol ~ 

■ Pu+l 



cfew(25+l-po^^) 
chw 

^^ shw{2s+2-po^^) 
sh2w 



Table 4. 



The symbol "hat" in table 4 stands for the Fourier transform. 

Before finishing this section, we want to point out that the (3 limit of the free energy ( |3.13|) 
provides the justification of the string hypothesis. Indeed, using results collected in table 2 and 
the formula ( |3.13|) , it is easy to obtain: 

lim/3y = -J2\dj{Q)\Hl + fJ-') = -IniV' (3.16) 

where N' = N {N/2) for odd (even). Equation above implies that the total number of states 
{N')^ is correctly reproduced by the use of the string hypothesis. 
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4. The Commensurable N-even Case 
N=2 SUSY 



From the results of the previous section (tables 3 and 4) we get the relevant information for 
unraveling the physics of the N even case. To make the resulting picture intuitively more clear, 
we will begin by considering the commensurable case, point A in figure 2. For generic N even 
this corresponds to the anisotropic antiferromagnet of spin j = and anisotropy 7 = ^• 



4.1 Commensurable Case: RSOS-structure 

A qualitative picture of the model can be obtained by comparing its string structure with that 
of the isotropic model with identical spin. In the isotropic case we have the same negative 
energy strings but a different pattern of the zero energy strings. In fact, all strings with length 
greater than 2j are allowed and all of them have zero energy. The elementary excitations for 
the isotropic model are holes in the Dirac sea of 2j-strings. The internal quantum numbers of 
these elementary excitations are of two types: a "vertex" -spin with value equal to 1/2 and an 
"RSOS"-spin defined by the Coxeter graph of the type A and Coxeter number 2j + 2. These 
two types of spin determine the dimension of the Hilbert space H^"-^ with n- holes [^: 

dimT^W = 2- Y: . . = 2-Y. — TT^(2cos(-— ))" (4.1) 

ii,...,i„_i=i i=i 3 + ^ ^3 + ^ 

with the incidence (2j + 1) x (2j + 1) matrix of the graph A2j+2'- 



K 



010 
101 








101 
010 



(4.2) 



The 2" contribution to ( |4.1| ) can be effectively obtained by ignoring all the zero energy strings 
with length smaller than 2j. In fact, the counting of the spin 1/2 degrees of freedom is. 
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essentially, identical to that of the XXX chain. The RS OS-contribution comes from the zero 
energy strings with length smaller than 2j. 

For the commensurable case (see table 4), we do not have zero energy strings with length 
greater than 2j. Instead, we have positive energy negative parity strings of length 1. By drawing 
analogy with an isotropic case, we may expect the same RSOS structure, but somewhat different 
"vertex" spin contribution. 

From the results collected in tables 3 and 4, we conclude (for even commensurable case) 
that the ground state is a Dirac sea of (2j = u — l)-strings and that there exist two types of the 
elementary excitations, degenerate in energy: holes of the Dirac sea and particles associated 
with the negative parity strings. The particle-hole structure of the commensurable chain is 
identical to that of XX model. These qualitative arguments indicate that the commensurable 
chain is equivalent to the free fermions equipped with RSOS-internal degrees of freedom. To 
prove this claim, we will proceed in step by step fashion: first, computing the dimension of 
the Hilbert space of n-elementary excitations, then, deriving scattering S- matrix and, finally, 
calculating the central extension. 

We start defining the spin Sz of the generic Bethe state \ip > as a difference in the number 
of roots with respect to the ground state 

S^d^jj >) = \^roots{\^lj >) - \^roots{\Gr. state >) (4.3) 

For even case we find with a help of the equation ( p.7| ) and the expression for the ground 
state density (table 4) the following: 

SM^>) = ^in.-ni%) (4.4) 
where n,^ is a number of positive energy strings with negative parity and n^i^li is a number of 



holes in the Dirac ground state. We can interpret ( [4.4|) as a renormalized spin and refer to 



-(h) 

^""^ as a spin of the state. Recalling that Po = and taking into account that is always 
an integer number, we obtain the restriction 

fly — n^uli = even integer {po odd) (4.5) 

Next, we proceed to compute the dimension of the Hilbert space with n-elementary excita- 
tions. The strategy will be to count all degenerate Bethe states containing the total number of 
particles (i/-strings) and holes equal n. From the hole-particle degeneration, at the commen- 
surable point, we readily get 2"-states. Additional degeneration will correspond to the zero 
energy strings with length smaller than 2j, i.e. the RSOS part. 

To count the RSOS degeneration we should, first, invert the equations ( p.7|) with a help of 
the inversion formulas presented in appendix B. For A^ even case we derive 

pi + Pi' = Si* py 
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Pj + Pj 
Pu~2 + Pul2 

pu-1 + Pu-1 



Si * {pfil + ; l<J<u-2 

Si * {pu-3 + Pu-1 + Py) 



(4.6) 



Pu + p. 



(h) 



Si * Pu-2 + 1 T 

Apol 



where Si was defined in section 3. 

Taking the Fourier transform of (^.61) and setting w = 0, we have for the non-negative integers 



nj[nj ') = pj{w = 0){pf'\w = 0)) the following system 



rii + n\ 

Uj + rij 
. ~{h) 



1 



ni% + k 



(4.7) 



with k = """^^^ ^ . Eqns. ( [4. 7] ) imply: 



(h) 



k = ^n,(z/-l-j) + ^^- 

j=i ^ 

v-2 



5]2n,(z-j)^^(z-j>0) + zni' 
i=i 



(4.^ 



From the eqns. ([4.8|) we get, for a given k, the set of pairs {ni,n-''^}, 1 < « < — 2, which 
characterize the RSOS degeneration: 



1/-2 



En 



i=l 



rii + n) 



Hi 



"count" 



(4.9) 



where 



2k. 



"count" = (ir )i,i = J2 



2i+i sin^f^^) 



i=l 



2i+2 ^ 
J + 1 



(2cos( 



^^2k 



2j + 2 



)) ; 2k = even integer (4.10) 



1 2j+i sin^ f ^' ) 

"count" = (K^w = E (-1)'^^ j+T ^^ '°'^2jT2^^'^ ' = '''^'^ ' ^ " 

(4.11) 
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In the formula ( [4.9| ) the sum is taken over the whole set of solutions to ( [4.8|) . 

Even though we don't have an analytical proof of the identity (|4.9| , [4.10| , [4.11D , we ran 
extensive numerical check on it, so that we are convinced in its validity beyond reasonable 
doubt. 

For 77,-elementary excitations {n = 2k) we get 

dim7^(") = 2" X "count" (4.12) 

Next, we will consider two possible cases: po-even integer and po-odd integer. The po-odd 
integer case corresponds to the spin j integer and the " count" number ( [4.101 ) coincides with the 
multiplicity of the identity in the tensor product of 2k copies of the spin 1/2 irrep, provided we 
use the truncated decomposition rule: 

min(ji+j2,2i-ii-j2) 

Jixj2= E 2j = z/-l (4.13) 

This is a well known decomposition rule for irreps of SU{2) Kac- Moody algebra with level 2j. 
The RSOS interpretation is obtained by associating with each state in H^'^^ a path of length 
n + 1 in the Bratteli diagram defined by the graph ^2^+2, starting and finishing at the zero 
point (see fig. 3). 

It is easy to conclude from fig. 3 that if n is odd, then the RSOS contribution is zero. This 
fact reflects the restriction (nj, — n^^\) = even integer we have mentioned above. 

The po even integer case corresponds to the spin j-half integer. In this case eqn. ( [4. 4)) does 
not impose any restrictions on the number of particles or holes of the state. For the formula 
(|4.11|) we obtain the following Bratteli diagram interpretation: 

i) 2k = even integer: Path in ^2^+2 start and finish at zero. 

ii) 2k = odd integer, k > j: Path in ^2^+2 start at zero and finish at j. 



4.2 S*- matrices 



The physics of the model will result from solving the following "inverse problem", known as 
higher level Bethe ansatz p2|. Given an arbitrary Bethe state \ip >, we should look for a 
scattering S-matrix such that the Bethe equations are equivalent to: 



e^'^^^l^ >= S{k,, A;,+i) . . . Sik„ %) . . . S{k„ kj_,M > (4-14) 

The S'-matrices in ( [4.14 ) will describe the scattering between the elementary excitations of the 
model. From the previous counting, we expect that the S'-matrix in ( 4.14 ) factorizes into the 
two pieces, each one associated with one of the two different types of spin. 
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The string structure of our model at the commensurable point (i.e. presence of only one 
string with index greater than 2j and particle-hole degeneracy) can be taken as indication 
that the spin 1/2 part of the S-matrix must be the one of the XX model. The following is 
a simple argument to make this claim plausible. In the isotropic case, the singlet and the 
triplet states differ by the existence of an extra zero energy string of length 2j + 1. This 
extra string is the reason for the difference between the singlet and triplet ^-matrices. In the 
commensurable case, we do not have this extra string and the "singlet" (particle- hole) and the 
"triplet" (particle(hole)-particle(hole), hole-particle) are practically the same from the point of 
view of the string content. 

To get the S-matrices, we will start with the system of eqns ( [4.6|) from which we can obtain: 



\^u~~i\ . , 4h) sh{u--2)vj . 4h) ^ , sh{v-2)w 



. J Pu-i + Pu-i - TTT u T^'^^'^ + ^^-i) + ~u T^P'^-s + • • • (4.15) 

ApqI 2chwsh{h' — l)w sh^v — l)w 

Note that the equation above has the same structure as the original Bethe ansatz equations 
(|3.71) , however, with one important distinction: the roles of and p^/*^ are completely reversed. 
Let us consider the simplest case: the Bethe state with two holes and one {v — 2)-string with 
rapidities Mi, U2 and "^"^"^ , respectively. The hole-hole S-matrix will contain two contributions: 
one, coming from the scattering between two (z/ — l)-holes and the other, from the scattering 
between (z/ — l)-hole and the (z/ — 2)-string. These contributions can be directly read from 



1 dw . sh{u-2)w 

\nSy-i^u-i{u) = -- / — smwu— — -— (4.16) 

2 J-oo w chwsh[u — \)w 



/CO 
dw cos wu 

The hole-hole S'-matrix is given by: 



shlu — 2)w 
oo shi^u — l)w 



S[u2i) = 5^-v-i(m2i)5^-i,i.-2(^) = e cH^sH^w ^ (417) 

Using the same technique, we can easily check that the particle-hole and the particle-particle 
S-matrices are identical to ( [4.17] ). This implies that the particle-hole system is, indeed, be- 
having as XX-chain with respect to the vertex spin. The formula ( [4.17|) gives us the RSOS- 
contribution to the scattering S'-matrix. For the case of the two elementary excitations the 
RSOS contribution is rather simple, in fact, it corresponds to the RSOS-Boltzman weight 

1 \ [ ] 
id i 
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for the A2j+2-niodel. 

Thus, we can represent our S'-matrices in the following symbolic form: 



Sxx S iisOS{p=l,r=2j+2) (4.18) 

Notice, that the S'-matrix we got, has the same structure as the S"- matrix for N = 2 integrable 
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S''=' = Sl='^Sl=l (4.19) 

where S^^"^ is the Sine-Gordon S-matrix for = 87r3/2 (which corresponds to g = 1 with 
q = e^'^/P and p = ^^^^2 ) • The N = 2 SUSY is not explicit in the Sine-Gordon model, but is 
realized through the quantum group symmetry in the nonlocal way [O, P|. 



4.3 Central Extensions 

Our analysis of the S'-matrices together with the counting of states led us to the conclusion 
that the commensurable chain is equivalent to the tensor product of an XX spin 1/2-chain and 
an RSOS model based on the graph A2j+2- This picture will be consistent if the central charge 
is given by: 

c = 1 + (2 — ) = — ^ ; 2j + 1 = - (4.20) 

^ 2j + 2^ j + 1 ' 2 ^ ' 

Two pieces which appear in the sum above are due to the XX and RSOS contributions, 

respectively. 

To see that this is, indeed, the case we now proceed to calculate the central extension from 
the low temperature expansion for the entropy. Following along the standard lines 0, we 
differentiate the system ( p. 15 ) with respect to u and compare the result with the system ( [4.6|) . 



Then, it is trivial to establish asymptotic validity of the following formula: 

Pj ~ — — — ; Vsound = Pol ; M > 1 (4.21) 

27[Vsound 1 + Vj 



Plugging (|4.21|) into the expression for the entropy S 



S 
L 



J2 J MiPj + pf^) Hpj + pV) - Pi Inpj - pf^ InpJ. 

i 

Y: J rfMp,(M)[(l + r/,)ln(l + r/,.)-r/,lnr/,] (4.22) 
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we obtain 



S^2TLf] - C{-^)} (4.23) 

where £ stands for the dUogarithmic Roger's function, defined as 

Taking ?7j,^j from table 1 and expfoiting "magic" formulas for the dilogarithmic functions ||33| : 

E.5L,/:(^) = if-2£(-i,) (4.25) 



£(x) + £(l-x) = ^ 



we, finally, obtain for 5* and the central extension c the following: 



^ 2'kTL ,\ 1 , , , 
4.26 



'^V sound dS _ 3'' 

ttTL df ~ ^ + 1 



c= ^^^^^ = ^ (4.27) 



Clearly, the formula ( [4.27| ) above is in perfect agreement with our prior anticipation ( [4.20| ). 



4.4 Integrable Deformations 

The XX-part of our model, by itself, admits nilpotent extension as we have discussed in section 
2. This extension corresponds to an addition of an extra magnetic field coupled only to spin 
1/2 degrees of freedom. The product of the nilpotent i?-matrix {N = 4) and an RSOS S'-matrix 

Rn'=4^ ^RSOS{p=l,r=2j+2) (4.28) 

describes our model away from the commensurable point (i.e. for generic spin). To check this 
claim, we, first, observe that the energies of holes e^-i and particles e^, presented in table 4 
are, essentially, equivalent to those of the nilpotent XX-chain. Second, the RSOS counting 
formula ( [4.9|) remains valid and so does our central charge result ([4.27|) . Note, however, that 
the particle- hole degeneracy is now irretrievably lost. 

We now explain the somewhat surprising result that the central extension does not show 
any dependence on our nilpotent magnetic field related to the irregular spin s. The point is that 
this internal magnetic field is selective in nature: it affects only fermionic degrees of freedom, 
but not RSOS degrees of freedom. Since, the central charge for the free fermions (with or 
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without magnetic field) is equal to 1 and an RSOS part is immune to the internal magnetic 
field, we conclude that all conformal degrees of freedom are preserved and, therefore, central 
charge does not change when we move away from the commensurable point. This is contrary 
to the behaviour of the higher spin XX Z model in the external magnetic field, which freezes 
out all parafermionic degrees of freedom and reduces the central charge from Pi^ZW^-value to 
just 1. 

We conclude this section with the following observation. We started by introducing the 
model based on the nilpotent representation of quantum group. This nilpotency somehow 
managed to survive the "hell" of the renormalization (filling in Dirac sea, etc.) and reemerged 
(like Phoenix) in the modified form in the formula for the physical S-matrices, as expression 
(|428|) indicates. 
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5. Odd Case: Beyond Bootstrap 



For odd, the commensurable case corresponds to (7 = = ^^y^)- Results collected in 
tables 3 and 4 of section 3, clearly indicate that the string structure of the odd case is a 
bit more involved than that of the even case, and the departure from the properties of an 
isotropic model (7 = 0,j = ^y^) is also more dramatic. First of all, the ground state strings 
are of length 2je// with jeff = | and not, as it could be naively expected, of length 2j. This 
circumstance should be interpreted as an indication that for an odd case our model is in the 
strong anisotropic regime. The strings with length smaller than 2jejj are of the zero energy 
which indicates the presence of an RSOS structure based on the y42j^^^4.2-graph. The pattern of 
the positive energy strings is also new. For the A^ even case, we have only the negative parity 
strings degenerate in energy with holes of the Dirac sea. For the A^ odd case, we have two 
types of positive energy strings: one with length 2j + 3 and another, negative parity string of 
length 1, designated as (z/ + 1) and //-strings in tables 3 and 4. At the commensurable point, 
the {u + l)-string and the holes of the Dirac sea are degenerate in energy. In this section we will 
describe some aspects of the system defined by the holes of the Dirac sea and the two positive 
energy strings in terms of the RSOS model with p = 2 and the restriction parameter r = 6 



The emerging picture is drastically different from the one we have found in the A^ even case 
where the model was shown to be equivalent to the tensor product of a vertex XX model and 
an RSOS model based on the ^2^+2 graph. Notice that the (p = 2, r = 6) RSOS model which, 
for A^ odd, replaces the XX-part of the N even case, describes the spin 1, S'[/(2)g-invariant 



Fateev-Zamolodchikov chain with = 1 ||34|| . 



5.1 Counting of States 

The first thing we would like to do is to count the degeneration of the Bethe states following, 
essentially, the same procedure we used for the N even case. Once again we invert eqns. (|3.7|) 
with a help of the identities presented in appendix B to get 
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- (S2 + Si) *pW + /(«)] ; l<j<u-l (5.1) 
(h) (h\ {~^Y COS 2770; ^ „ ^ 

Here, notations are the same as in section 3. Applying the Fourier transform to the system 
above and setting w = 0, we derive, after a bit of labor, the following: 



ni + n\ = 


2 ^ 


rij + ny = 









l<j<u-2 (5.2) 



where 



'k = hy ^ — ; 2k = even integer (5.3) 



Comparing system (5^) with system ( [4.7| ), one cannot help noticing that they are remarkably 
similar. In fact, the relation between them is simply 

k ^ k (5.4) 
z/ z/ + 1 

Note, however, that unlike 2k, 2k is always an even integer number. 

Repeating the same analysis as in section 4, it is straightforward to infer that the total 
number of the degenerate Bethe states with the same {k, fij,} is a product of two factors. The 
first RSOS factor is due the to zero energy strings contribution and is equal to the number 
of paths of length 2k in the y42j^^^+2-Bratteli diagram, which start and end at zero point. 

The second factor 2^(^^""') is related to the degeneracy between holes and (z/ + l)-strings. 
Note that this degeneracy disappears once we move away from the commensurable point and, 
consequently, the second factor reduces to 1. 

Again, we can describe RSOS factor as the multiplicity of the identity in the tensor product of 
2k-copies of the 1/2-irrep. This tensor product is defined by the truncated decomposition rule: 

min(ji +j2 ,2je// - ji -32) 

31 X J2 = J2 3 (5.5) 

The explicit formula for the RSOS factor (k-fixed) is the same as ( [4.10| ) with je// substituted 
for j and k is substituted for k. 
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For the spin of the Bethe state >= \n^, n^+i, rL-yl2 > one finds 

Sz{\^ >) = '^rootsd^ >) - \\roots{\Gr. state >) = Y(ri,.+i - (5-6) 

This imphes (for odd) the restriction 

(fii^+i — = even integer (5.7) 

From ( |5.6D we also observe that the negative parity strings do not contribute to the total spin 
of the state, which indicates that the //-string behaves as a spin "singlet". 

Ignoring, for the time being, the zero energy strings, we now proceed to map the Bethe 
states |n,^, n,^+i, n^^2 > o^i^o the states of the (p = 2,r = 6) RSOS model. This model is 
defined by the disconnected graph (see fig. 4). We will limit ourselves to the states defined 
by the first graph I. In terms of this graph we can identify the holes and u + 1-strings with 
the 2-point links (ji,j2); \ji — j2\ = 2 and the //-string with the closed link (3,3). We will 
also impose on the RSOS states periodic boundary conditions. It is easy to check that this 
identification is compatible with the counting and restriction (|5.?1 ) we have described above. 
In fact, for the given values nu,nu+i and n^J^l2^ the corresponding number of the RSOS states 
related to graph I is 2'-"'"+^ "+2''. In this way we capture the "singlet" nature of the //-particle. 
Restriction (|5.7|) can be easily translated into the RSOS language as a requirement that one 
should consider only RSOS states which start and finish at the same point on the graph. If 
we use graph I to describe holes and positive energy strings then we will find, as a basis of the 
Hilbert space with the total number of n-elementary excitations, the following: 



where 



{31,32, ■ ■ ■ ,jn+i) <^iJi, Jn+i) (5.^ 



n = n^ + n^+i + 11^1^12 3i=3n+i jG {1,3,5} (5.9) 

N = 2hy + fiyj^i + n^^l2 Ji = Jn+i = J ^ {^, ^, ■ ■ ■ , Jeff} 



and 



2 for {n,y^i + ri^+2) couples 
and ji = 3 for couples 



■1 



States ( |5.8| ) for the different partitions of n into the sum of the three positive integers provide 
a basis of the Hilbert space with n-elementary excitations. 
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The basis (|5.8| ), which is aheady compatible with our counting, imphes that the S'-matrix 
for the elementary excitations is a product of the two RSOS S'-matrices. In the next section, 
we will show that this is, in fact, the case. Before passing to this issue, let us briefly comment 
on graph II. This graph can be interpreted as describing strings in a discrete target space 
characterized by the graph ^3 According to this graph, we can suggest the following 

"stringy" interpretation. The holes and {u + l)-strings can be associated with "closed" RSOS 
configurations (2, 2) and (4, 4), and the negative parity string u with "open" RSOS configuration 
(2,4). If we maintain this "stringy" interpretation of the graph II, holes and {u + l)-strings 
will appear as bound states of //-strings. 



5.2 S-matrices and Bootstrap 

Let us consider the (p = 2, r = 6) RSOS model which is a descendant of the (p = 1, r = 6) model 



and is obtained by the standard fusion procedure This model can live in two different 

regimes: 

regime (I): The ground state is a Dirac sea of length p = 2-strings and 
there are no positive energy strings. The central extension 
is given by: 

3p 2(p + 2) . . . .^ 

p + 2 r[r — p) 

regime (II): The ground state is a Dirac sea of length r — 2 = 4-strings. 

Strings with length 1, 2, 3 = (r — 3) are of the positive energy. 
Holes are not allowed in the (r — 2) Dirac sea. 
The central extension is given by: 

c = 2- - = l; r = 6 
r 

The elementary excitations in regime (I) are holes of the p = 2 Dirac sea, while in regime (II) 
they are associated with the positive energy strings. Notice that for (p = 2,r = 6) the central 
extension is the same in both regimes. To match the dynamics of (z/ + 1, i/)-strings and holes 
of our model we will choose regime (II). 

The Bethe equations of the (p = 2,r = 6) model can be written as follows |23|: 



pf'+i:A^i*pk-A<j<^ (5.11) 



k=l 



with the functions aj^l and A^!^l best defined by their Fourier transforms: 
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2ch\ 



^iS(^) = '2aj{w)nk{w) , j>k 



(5.12) 



ttj [W ) 



sh{2w) 



nfc(iy) = ct/i(-)s/i(— ) 



To compute A^^^l for j < k we should use the symmetries of the model 



^(4) _ ^(4) 
^1,2 — ^3,2 

= m (5-13) 

^^(4) _ ^(4) 
^1,3 — ^3,1 

In the regime (II), the positive energies of 1, 2, 3-strings are given by: 

e, = (5.14) 
Using (|5.12|) and ( ^.13[) we conclude that the 1 and 3 strings are degenerate in energy 



cft(f) 

€2 - 



chw 

This will agree with our counting if we identify the 1 and 3 strings of the (p = 2, r = 6) RSOS 
model with the holes and u + 1-strings. However, the energies do not match in a peculiar way. 
From table 4, we get (up to the overall factor) 

eu+2 = Gu+i = , (5.16) 
chw 

1 

e„ -- 



chw 



These are the values for the commensurable point for arbitrary odd. Comparing ( |5.15| ) and 
(|5.16|) we, first, observe that for the (p = 2, r = 6) RSOS model the string 2 may be interpreted 



as a bound states of two 1-strings. The u + 2-holes and u + 1-strings have the same energy 
as 2-string of an RSOS model, however, the negative parity //-string has twice the energy of 
the 1, 3 strings. Furthermore, it is straightforward to show that this mismatch of energies 
persists through the whole hermiticity interval ( |2.35|) . Before we give an interpretation of these 



results, let us see how S'-matrices, or better to say, the scattering shifts of the odd case and 
(p = 2, r = 6) RSOS model fit together. 
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The scattering shifts for the interchange of j and fc-strings are related to the kernels Aj^l in 
(|5.11|) . The (p = 2, r = 6) RSOS model is a descendant of the (p = 1, r = 6) RSOS model which 
completely determines the Bethe strings scattering. The explicit dependence on the value of 
p = 2 only appears in the energies ( ^.151) i.e. in the left-hand side of the Bethe eqns 
Using (|5l2| , |5ll we find 



Atl - 1 = A^'^ 



3,3 



^(4) 
^2,2 



1 — A^^^ — 4^'*-' 
— ^1,3 — ^3,1 

1 



(4) 



1 



1 + 



chw 



^2,3 — ^2,1 — ^3,2 — ^1,2 



2chw 



chw 



(a) 
(b) 



(5.17) 



It is important to notice the existence of the extra symmetry reflected in eqns. (a, c) above 
which is peculiar for the (p = 2, r = 6) model. Using the results above it is trivial to obtain for 
the S'-matrices the following: 



(5.1^ 



5i(3),i(3) = tan(f + zf) 
^2,2 = tan\l + t^) 

Rewriting the Bethe ansatz eqns. ( ^.71 ) in such a way that roles of p,y+2 and p^^l2 are reversed, 
we get 



4poI 



+ 2ch{ 



chw 

1 s/i(z/ — l)w 

chw shvw 
w sh{u — l)w 



sh(u — l)w., jh) 



shiJW 
sh{u — l)w 
chw shvw 



shvw 



, ,w. shiv — 2)w „ 

+ 2ch{-)^- '—p,_2 + ... 5.19 

2 shvw 



Pu+2 + Pu+2 + ( 



+ 



ch{ 



)(pS 

sh{v — l)w sh{v — l)w 



2 chw 



1 - 



sh{v — l)w, 
2chwshvw 

-)pu + 



chw shvw shvw 

Reading off the scattering shifts from the r.h.s. of (|5.19| ), we have 



Pu+l) + 
pu-1 + • • • 



SholeMole{u21j 



SuAu2l) 



X 



dw ■ / .sh(iy-i)w ohJL(li2L 
— sm(ii2i«')— — r-^ 2;y V 2 

w V / chwshuw ^ 



shf^\ 



— X tan{ — h i 

^ + i 4 4 ■ 



p ■ / \ sh(iy—l)v 
■ SimUoilVj—rr TT^ 



(5.20) 



X tan^i— + i 

^ A 



IT 



^hU"^ - shf^{^ - 10 shl^{^ - |z) 



TIU21- 



'2u' 
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Here, the first expression describes an evolution of tlie state consisting of two {u + 2)-holes 
and one (z/ — l)-string with the rapidities ui,U2 and "^^"^ , respectively. The second one 
corresponds to the state which has two //-particles with rapidities (mi,U2) and one {u — 2)- 
string with the rapidity "^^"^ . Let us summarize the salient features of the equations above. 
They, clearly, indicate that our S-matrix is a product of two distinctly different S'-matrices. 
The first one is related to the RSOS defined by graph A2j^jj+2- The second one is due to 
(p = 2,r = 6) RSOS model. Note that the second RSOS comes without the usual "vacuum- 
polarization" contribution. Also note that //-string behaves as a bound state of two holes. This 
is not consistent with our previous result ( |5.16[ ), thus signaling the breakdown of the bootstrap 
procedure. 

Once we have shown the matching of the scattering shifts between the odd model and 
the product of the two RSOS models: (p = l,r = 2je// + 2) = 2, r = 6) we can return 
to the question of the energies mismatch we have left open. The reader should be aware that 
the Bethe study of the scattering shifts we are performing here, is not sufficient to uniquely 
determine the dynamics of the holes, v + 1 and v strings. In fact, the scattering shifts will be 
the same for (p = l,r = 6) and for (p = 2,r = 6) RSOS models. Our strategy is to get pieces 
of the puzzle. First, we observe that the RSOS model with the restriction parameter r = 6 
correctly captures the scattering shifts. Second, use of (p = 2, r = 6) RSOS- descendant model 
leads to the correct counting of the Bethe states, gives the right central extension, however, 
fails to reproduce energies of the three types of strings which define the spectrum of the N odd 
case. Two comments are now in order. First of all, a mismatch between bootstrap properties 
of energies and of scattering shifts is typical of the nonfusion models and, therefore, we can 
imagine that the model we are seeking is some nonfusion descendant of the (p = 2,r = 6) 
RSOS model. Secondly and perhaps more likely, we can attempt to make contact with the 



Z3 parafermionic extension of the N = 2 SUSY S'-matrices defined in [16|, i.e. to match our 



spectrum with the one of the affine Toda SU{3) (also see comment at the end of this section). 
If this is the case, we arrive at the nice physical picture that the commensurable spin chains 
for the even and A^ odd cases provide a useful framework where to study the Z{2) and Z{3) 
N = 2 integrable models. Finally, we would like to draw the following important comparison. 
The most fundamental difference between A^ even and A^ odd cases is that for the A^ even, the 
commensurable point sits just on the border of the weak anisotropic regime, while for A^ odd, 
the commensurable point lies in the strong anisotropic regime. The analogous difference exists 
between the massive N = 2 model perturbed by the least relevant operator and N = 2 model 
perturbed by the most relevant operator. To complete this section we move to the computation 
of the central extension. 

Taking advantage of the relevant results from table 2, formulas ( [4.23| , ^.25| ) and additional 
identity 
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we have for entropy S 



^(7rT^) = 2^(^)-2A^) (5.21) 



T^V sound ^2 ^ «^ (1 



and for the central charge 

The RSOS representation of our model 

RSOSp=2,r=6 RSOSp=i^r=2jeff+2 (5.24) 

gives us the central extension 

^ = (2 - |) + (2 - ^-^) = , (5.25) 

6 2je// + 2 Je// + 1 

provided that both RSOS models are in the regime (II). 

Clearly, the formula above is in perfect agreement with our exact result (|5.23|) . 

We would like to finish this section in somewhat speculative vein to express what may 
be the most natural extension of even case results to odd case (or, more generally, to 
generic g-root of unity case). Since even case describes free fermions "dressed" by RSOS 
degrees of freedom, then it is plausible that the generic g-root of unity case corresponds to free 
parafermions upgraded by RSOS degrees of freedom. These free parafermions are "free" in a 
sense of possessing additional symmetry related to the Onsager algebra (see ref. |^ for the 
latest review). We hope to expand on these matters in our future publications. 
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6. Magnetic Properties of Nilpotent 
Spin Chain 



In this section, we briefly discuss behaviour of our nilpotent spin chains in the external magnetic 
field [h). In ref. p2[ we have shown that at h = T = our system exhibits phenomenon of 



ferrimagnetism, that is, the ground state spin is greater than zero-antiferromagnetic value, but 
smaller than maximum possible ferromagnetic value. From the discussion given in the previous 
sections, it is abundantly clear that this nonvanishing magnetization is due to the presence of 
"internal magnetic field", which selectively effects available degrees of freedom. By placing the 
system in the external magnetic field, one may hope to demonstrate intricate interplay between 
these two fields, and to "probe" the model beyond extreme long distance universality, which is, 
to extend, nilpotency "blind". 

Let us choose "up" small magnetic field in such a way that the ground state Fermi band 
shrinks from — oo < m < oo to —B < u < B where -B ^ 1 is the Fermi momentum. Following 
Yang and Yang , we have for T = 



roo roo 

PjJu + B) = p^{u) + J{u - fi)p+{p)dfi+ J{u + p + 2B)p^{p)dp (6.1) 

JO JO 

where is an unperturbed density describing the distribution of the ground state strings (see 
table 4), J{w) = — \ . — 1 and p+{u) is density of holes induced by the external magnetic 
field, and which are located outside of the completely filled Fermi band. 

In the eqn. ( |6.1| ), the last term appearing in the right hand side can be treated as a 
perturbation, provided B ^ 1. Iterating (|6.1|) , it is straightforward to obtain 

p4«) = p^t\n)+P^l\u) + ... 
p,,{u + B) = pf («) + /- J(«-p)pf(0)(/.)ci/. (6.2) 
. - /o°° J(m + p + 2B)pf^ (p)dp = pj ^ (u) + /o°° J{u- p)pV^^ (p)cip , etc. 

and, for the magnetization M, 

M = s- ^{Bethe roots) = Mq - 2(1 + J{0)){pf\o) + pJ^(O) + . . .) (6.3) 
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where, as usual, "hat" designates Fourier transform and Mq = s — Pj^iO), which represents the 
nonvanishing ground state magnetization {h = 0), was given in ref. [^. 

One can recognize (|6.2|) as Wiener- Hopf eqns. which can be solved in standard fashion to give 



< Wn(j)=l < Wn{j)=2 < ■■■ 

Here, a„ = reSu,=iwnPjo{w)] bj = reSw=iu,j+{w) and d.p. stands for the double-poles contribu- 
tion which is not important for the following. Kernels J±{w) are defined as 



J^{w) = Jz\w)- l + JH = itM; j^(oo) = l (6.5) 
Explicit formulas for J+ are given below: 

j^(^) = yigg.f(poinpo-(Po-i)in(po-i)) X ; Neven 

i+H = yge^fboinpo-(po-|)in(po-i)-inV2) ^ r(^O^jgf) . N odd {N > ^) 

(6.6) 

Note that J±(w) has singularities only in the upper (lower)-half plane and is analytic and 
non-zero in the lower (upper)-half plane. 

Before we move on, two comments are in order. 
First, the Fourier image Pjo(w) has two kinds of poles: 

a) those, which are present for s = j (regular spin) and 

b) additional ones, which appear only when s ^ j (i.e. in general nilpotent case). 

It is precisely through the "b"-poles that nilpotency manifests itself in the magnetic character- 
istic. As will be demonstrated shortly, the "b" -contribution is not invariant under the change 
of sign of the external magnetic field. 

Second, this comment concerns the pole structure of J-t (w)-kernels. The reader can easily 
convince oneself that this structure depends on whether u is even or odd. Even though the 
difference in the analytical properties is not that important for small h, it may become the 
relevant factor, once the strong magnetic field is imposed. Note in this regard, that the Hilbert 
space of the excitations (for even) is, indeed, sensitive to the //-parity, as was shown in section 
4. 
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To proceed further, one must establish the relation between Fermi momentum B and the 
magnetic field h. The easiest way to do this is to impose a condition that the energy of a hole, 
located at the Fermi surface, vanishes, i.e. 

e_^(n = 0) = (6.7) 

where is the energy of a hole, subject to 



— ApqI pjJu + B) + zh = e^{u) + / J{u — fi)e^{fi)dfi + / J{u + + 2B)e^{fi)dfi (6.8) 

JO Jo 

with z = po(^) for odd (even). 

Eqn. ( |6.8| ) above is very similar to the eqn ( |6.1| ) and can be iterated in the same fashion to 
yield 



e+{w) = -4po/(pV^(w) +pW) + •••) + T7 r^TT^ + ^ -^TTT^^^"^^ + • • ■ 

(6.9) 

In terms of €^{w) one can rewrite ( |6.7|) as follows: 

lim iwe+{w) = (6.10) 



Ji)— >oo 



It is reasonable to expect that the effect of the interplay between "internal" and external 
magnetic fields will be visible in the lowest order of the perturbation theory, if the number of 
"affected" degrees of freedom is comparable with the total number of degrees of freedom, i.e. 
if po is not too large. To show that this is, indeed, the case let us calculate magnetization M 
for po = 5. 

Making use of (|673| , |6^ , |6.1CI|) one finds 



h ,/i,2 2 , r2(l) 32sin(7r5s), ,^ ^ ^ , 

M{h) = Mo + — : + {jf-rr-2 rT-\^ v(^\vi 5x + ^ -) ■ ■ ■ 6.11 

27r/ / 457r^ cos(7rds) r(|)r(— |) vr 

where 5s = s — j . 

In the discussion above it was assumed that h > 0. The best way to treat h < case is to 
exploit "charge conjugation" authomorphic properties of our nilpotent hamiltonian H, namely 

H{6s,-h) H{-6s,h) (6.12) 
s N' -1-s 
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Note that these properties are the direct consequences of ( p.34| ). 
Now it is obvious that 



M{-\h\)=Mo + ^ + (}^ 



457r2cos(7r(5s) 'r(f)r(-f; 



32 sin(7r5s) 



TT 



} 



(6.13) 



Comparing eqns. (|6.11|) and (|6.13|) , one cannot help noticing that the coefficient of the quadratic 



term does depend on the sign of the magnetic field. Clearly, this effect bears the signature of 
the irregular irrep., since it only exists when the systems is away from the commensurable 
point. For po ^ 1 one should keep sufficiently many terms in the perturbation expansion to 
"catch" this highly nontrivial effect. 

Finally, let us present M{h) for the generic Pq ^ 1 

M-Mo = —- {1 + i ; ' ) -r..^ g} + ... (6.14) 

2poIWn=l AipQlJ_{lWn=l)an=l 

where g = '"-^f-^^"r'^''t • 

It is interesting to observe that the ffist term in the expansion ( |6.14[ ) is always universal (which 



is consistent with universality of the central charge) and that next to leading term depends 
on 5s. The presence of the fractional powers of h in the formula above should be taken as an 
indication that M has nontrivial analytical properties. 

To conclude this section, we would like to comment that the usual criteria of the small 
magnetic field h <^ I does not quite work in our case. Indeed, note that h in the second term 
of ( |6.14| ) appears in combination j^-z^^jj^- Keeping h fixed and varying s towards the orbifold 



point can make this combination very large. In fact, it can be shown that for arbitrary small 
J, there is a region adjacent to the orbifold point where the system is in completely ordered 
ferromagnetic state. We hope to report on the strong magnetic field properties of our system 
in the future publications. 
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Concluding Remarks 



In this paper we have analysed the special features of integrable higher spin XXZ chains 
resulting from imposing spin-anisotropy commensurability. This study can be considered from 
different points of view each one providing a direction of research we feel worth pursuing. We 
sketch some of them currently under study: 

i) Generalization of the concept of free fermions and N = 2 SUSY. Starting with the work 
of Onsager free fermions are the basic tool to solve a large class of 2D statistical models. The 
simplest generalization is to work with GI{1, l)g-invariant chains or equivalently with nilpotent 
deformations of S'f/(2)g-chains with = 1. We found for commensurable chains with = 1 
{N even) that the spectrum of excitations consists of free fermions, "dressed" by the RSOS 
degrees of freedom, with their "vertex dynamics" completely determined by the SU{2)q {q^ = 1) 
quantum symmetry. This symmetry underlies a new non-local realization of = 2 SUSY. The 
connection between free fermions and N = 2 models may shed some new light on the Painlevee 
structure recently found in the study of = 2 models ||39[. Our results for the = 1 {N odd) 
chains show the way to find a definition of free parafermions. Some hints of the non trivial 
structure of such systems already appear in the observed failure of naive bootstrap. 

ii) Integrability in higher dimensions. N = 2 SUSY theories, free fermion models and their 
generalizations are natural candidates for a direct 3-dimensional interpretation. A preliminary 
research in this direction was started in |^D|. The integrable deformations associated with 
nilpotent irreps we have considered in this paper allow, in some simple cases, the interpretation 
in terms of dimensionally reduced 3D-models. 

iii) Selective magnetism. Our models can be interpreted in terms of an "internal" magnetic 
field which selectively affects available degrees of freedom. As a consequence, the response of 
these systems to the external magnetic field exhibits asymmetry under the change of sign of 
the external field. Although, this effect is relatively small, it can be taken as a signature of an 
extra dimension in the spirit of comment (ii). 

We believe that the study of the commensurable spin chains provides good opportunities to 
link a variety of topics which are ripe for unification. 
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Appendix A 

Hermiticity Regions for Nilpotent 
Hamiltonians 



In this appendix we shall describe the intervals of the generalized spin s for which the corre- 
sponding hamiltonian with anisotropy 7 = 2ttM/N is hermitian. We shall assume that and 
M are relative primes such that 1 < M < A^. Then the problem is to find the regions where 
the variable x = 2s {0 < x < N) satisfies 

M M 
sin2n — (r + 1) sin27i — (x — r)>0 (A.l) 
N N 

for r = 0,1,..., N' -2. 

Considering a few examples, one extracts the following facts: 

i) There are M intervals with positive parity and M intervals with negative parity which 
alternate. 

ii) The width of each interval {xi,X2) is given by: 

= (A.2) 

^ 2/M N even ^ ^ 

iii) The orbifold points Ai = g^^ and A2 = q^^ that delimitate the intervals are related by 
charge conjugation: 



A2 = (Ai) 



DC 



(A.3) 



From these data we want to find the values of Xi and X2 which, in turn, imply the patterns 
of breaking of the nilpotent irreps at the orbifold points. From ([A.2|) and (|A.3|) we get: 



Nodd: kN-{l + Xi)2M =1 
Neven: |A^-(1 + xi)M =1 (A.4) 
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where k is some integer. These equations admit several solutions for xi which differ by Axi = 
We shall choose the solution for which d = 1 + xi gives us the dimension of one of the 



irreps at the orbifold point. The other is equal to A^' — d. Both cases in ( |A.4D can be treated 
simultaneously if one recalls the definition of A^' and one defines M' as: 



N odd : N' = N, M' = 2M' 
N even : A^' = f , M = M' (A.5) 

Using (|A.5|) , equation (|A.4|) becomes: 



kN' -dM' = 1 (A.6) 
Since A^' and M' are relative primes, this equation has a unique solution for k and d (both 

N' _ n. 
M' 7" 



integers) which can be found using the fraction expansion of po = ^ = ~ 



A^' 1 
Po = jp = ^i + T77^ (A.7) 

The integers z/j(i = 1, . . . , a) can be computed using the Euclid algorithm. In terms of them 
one defines the reduced fractions: 

1 1 A^' 

^1 = T^i, 52 = vi^ , 5■i = Vl^ . . . , 5„ = — (A. 8) 

which can also be written as 

5^ = ^- (A.9) 



where iji and computed from the recursion formulas: 

Vo = 1, 2/1 = ^1, Vi = J^iVi-i + yi-2 

Zo = 0, Zi = Z/i, Zi = ViZi^i + Zi_2 (A. 10) 

Notice that Ua = N' and z^ = M'. One can easily prove by induction: 

yiZi-i - ZiVi^i = {-ly ,i>l (A.ll) 
Setting i = a in ( |A.11| ) we get: 



N' z^^^ - M' y^., = {-ir (A.12) 
and comparing with ( [A.6|) we, finally, obtain: 
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a even : k = Za-i, d = Va-i 
a odd : k = M' — Za-i, d = N' — i/a-i 



(A.13) 



Therefore, for a even or odd, the dimensions of the regular irreps into which the nilpotent 
irrep breaks are i/a and A^' — i/a- Furthermore, we shall now show that both numbers i/a and 
A^' — i/a are Takahashi numbers. Indeed, the Takahashi numbers are the pairs {nj,Vj) which 
give the string length and parity of the allowed strings for a given anisotropy 7 and can be 
computed using the formula [[1^ : 



-J - Vi-i + U - "n^dVi for mi<j< rrii+i 
where mo = 0, mj = X]fc=i and the parities Vi are: 



(A. 14) 



1, = -1, Vj = exp in 



Po 



for j ^ l,mi 



(A. 15) 



Here [x] means the integer part of x. For rational values of po the index i in ( |A.14| ) runs 
from to a and one has an an extra condition besides ( [A.14| ): 



rij < i \fj (A. 16) 



The allowed values of rij given by ( |A.14| ) form the i-th Takahashi zone. Thus, there are 



a + 1 zones. Condition ( [A.16| ) implies that the last zone, i.e. i = a, has only one string of 
length = Va-i- On the other hand, using (|A.10|) one gets for the last string of the zone 
i = a - 1 the length n^^.i = A^' - y^-i- 

We have proved that the lengths of the two last strings in the Takahashi-Suzuki labeling 
coincide with the dimensions of the regular irreps into which the nilpotent irrep breaks at 
the orbifold points. This situation is quite different from the one considered by Kirillov and 



Reshetikhin |p (see also ref. |26|) in their study of higher spin XXZ models. They study 
situations where the anisotropy 7 is such that 2j + 1 is a Takahashi number, where j stands 
for the spin of the chain. We observe that in our case 2j + 1 = N' is not a Takahashi number, 
although, the irreducible pieces into which it breaks at the orbifold points are, indeed, the 
Takahashi numbers. 

Below we show some patterns of breaking for low values of A^. 
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A' 
4 
"3" 
6 



,1/ 
1 

"1^2" 

1, 5 



1,3,5,7 



Breaking 



1 e 1 

2 e 1 
2 e 1 



5 
5 
10 
10 



1,4 
2,3 
1,9 
3,7 



3 2 

4 e 1 
4 e 1 

3 2 



12 
12 



1,11 
5,7 



5 
5 



Table 5. 
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Appendix B 



The Inversion Formulas for Kernels Ta 



Starting with the formulas ( |3.8|) , it is straightforward (though tedious) to verify the following 
inversion relations exhibited in the table below. 



even 


N odd 


- 2chwi^J+^,k + ^j-l,k} + ^j,k 2chw ' 

i<j<jy-'2. 

-^u,k — 2chw '^~^^k + ^u,k 


^j,k - 2cL{"^i+i,fc + 4-i,fc} + Sj^k ; 1 < j < 2 

^u,k — 2ch^ l^^u-l,k ^^u+L,ki 1 ^u,k 2ch^ 
-^v+l.k — 2ch^ '^^k 5v+l,k 
Ay-\.2,k — 2ch^ '^'k + '5(/+2,fc 



where u 



Table 6. 

N /■N-i\ gygj^ (^j\f odd); symbol "hat" above kernel A designates Fourier 



2 V 2 

transform and, finally, the connection between kernels A and T is given by: 



Aj^kiu) = Tj^kiu) + {-lY'Sj,kS{u) 
The inversion formulas presented in the table above are quite similar to those discussed 



in ref. [0. The main difference is related to the fact that ref. deals with the generic 
anisotropy 7 and as a result infinitely many Takahashi zones enter into the analysis. Since in 
our case, g is a root of unity, the number of relevant Takahashi zones is always finite. This fact 
is reflected in the presence of extra 6 function pieces which appear in our inversion relations. 
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Figure Captions 



Fig. 1. Hermiticity Regions 

Fig. 2. Spin-Anisotropy Plane and the Commensurable Point 

Fig. 3. Bratteli Diagrams for ^2^+2 graph and the Paths for po even (odd) 

Fig. 4. (p = 2, r = 6) RSOS graph 
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